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' Abstract. We present high- precision Monte Carlo data for the phase diagram of 

^ . a two-species driven diffusive system, reminiscent of traffic across a narrow bridge. 

Earher studies reported two phases with broken symmetry; the existence of one of these 
has been the subject of some debate. We show that the disputed phase disappears for 
, sufficiently large systems and/or sufficiently low bulk mobility. 
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29 ' 1- Introduction 

o 

^ . 

O ■ When driven far from thermal equihbrium, interacting many-particle systems display 
. a wide variety of unexpected behaviours [H Ej. Generic long-range correlations, new 
phase transitions and universality classes, and anomalously fast coarsening are amongst 
. the most striking phenomena. In contrast to equilibrium systems, our notions of 

Ch , energy-entropy competition fail so that, even in low spatial dimensions, a plethora 
C_> , of novel phases and phase transitions have been discovered 12 El- In the absence 
of a fundamental theoretical framework for nonequilibrium systems, these transitions 
^ '. are usually identified with the help of Monte Carlo simulations, and subtle finite-size 
analyses are essential in order to establish their nature in the thermodynamic limit. 
Only in a few cases can such questions be settled by recourse to exact solutions. 

Asymmetric exclusion processes (ASEPs) constitute a particularly simple class of 
far-from-equilibrium model systems 0. Over the past two decades, they have emerged 
as a kind of "laboratory" in which nonequilibrium phase transitions can be studied. 
ASEPs have also been invoked to model a variety of physical systems, such as biological 
transport, traffic problems and gel electrophoresis. Defined on a regular one- dimensional 
lattice of length L, particles and holes perform nearest-neighbour exchanges biased in a 
specific direction. If exchanges in the suppressed direction are completely prohibited, the 
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process becomes totally asymmetric (TASEP). For periodic boundary conditions, the 
TASEP stationary state is uniform; in contrast, if particles are inserted and extracted 
at the boundaries with certain rates (a and /?, respectively), three distinct phases are 
observed: a low-density (a < 1/2, a < /3), a high-density {a > (3, P < 1/2), and a 
maximal current (a > 1/2, (3 > 1/2) phase jHEllSl- 

If two species of particles are present, phase transitions can occur even in periodic 
systems [3 El- If the two species (labelled "positive" and "negative") are biased in 
opposite directions, move into vacant sites (holes) with unit rate, and exchange with 
the opposite species with another rate (7 < 1), then solid clusters of particles separated 
by nearly empty domains form easily. Though this particular model is proved to 
display only a homogeneous phase Pj, dramatically different properties emerge when 
it is generalized in a seemingly trivial manner - from 1 x L to 2 x L, modeling a "ttfo 
lane road" (see jTHI and references therein) . 

In this Letter, we consider the same biased two-species dynamics, but with open 
boundaries. In this "bridge" model, positive (negative) particles enter the system at 
the left (right) with rate a and exit at the right (left) with rate (3, reminiscent of traffic 
crossing a long road bridge in both directions. While cars move freely on clear stretches, 
they slow down for on-coming traffic, reflected in the bulk mobility 7. This dynamics is 
clearly symmetric with respect to charge-parity (CP) transformations. The flrst study 
[TT] of this model focused mostly on 7 = 1 (i.e., particle-particle and particle- hole 
exchanges occur with the same rate in the bulk) and reported four phases. Two of these 
phases are characterised by CP-symmetric steady-state particle densities and currents, 
while the other two spontaneously break the CP-symmetry. Speciflcally, choosing a = 1 
and lowering (3, the system undergoes a transition from a symmetric (the so-called 
"power-law") phase with maximal current, equal particle densities and no holes to a low- 
density symmetric (Id-S) phase where the particle densities are still equal but below 1/2. 
As f3 is lowered further, spontaneous symmetry breaking occurs: While both densities 
remain below 1/2, the two species are associated with different currents and densities. 
At even lower (3, this "low-density low-density" [Id/ld) phase undergoes a flnal transition 
into a "high-density low-density" [hd/ld) phase in which one of the particle densities 
exceeds 1/2. While the existence of the two symmetric phases is supported by an exact 
solution for (3 = 1, evidence for the two asymmetric phases is restricted to a mean-fleld 
theory and simulation data for fairly small systems. The existence of the (W/W)-phase 
was subsequently disputed jT2] , based predominantly on simulations of a suitably deflned 
free energy functional which can be extracted from particle density histograms. Arguing 
that this functional is not very sensitive to the subtle differences between the (Id/ld)- 
and {hd/ld)-ph.ases, more accurate data for particle density histograms were analysed in 
|13j . providing further evidence for the (/(i/W)-phase albeit with some unusual flnite-size 
properties. 

In this Letter, we investigate this model further, performing high-precision Monte 
Carlo (MC) simulations for much larger system sizes and an expanded range of 7, with 
< 7 < 1, in an attempt to clarify the nature of the (W//(i)-phase. We present very 
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strong numerical evidence that only a single symmetry breaking transition persists in 
the large L or small 7 limit, in the sense that the range in j3 over which the (Id/ld)- 
phase is observed narrows until it becomes unresolvable. As in the periodic case [211, the 
behaviour of low- dimensional systems can easily be misextrapolated unless very careful 
finite-size analyses are performed. Our findings confirm how subtle finite-size effects 
can be when studying driven diffusive systems. The importance of exploring larger and 
larger systems cannot be overemphasised. Since a reliable extrapolation of MC data is 
essential for the analysis of many nonequilibrium phenomena, the understanding of such 
peculiar finite-size effects may have repercussions well beyond these simple models. 

This Letter is organised as follows. We first introduce the model and our simulation 
method, then summarise our data, and conclude with a brief summary. More details 
will be reported in [T^ . 

2. Model and methods 

The model is defined on a one-dimensional lattice of length L. An occupation variable 
Si = is assigned to each site i = 1,...,L, indicating its state as empty, (0)i, or 
occupied by a positive, (+)«, or negative, (— )«, particle. Positive particles enter the 
system on the left with rate a and move to the right with rate 1 if their neighbour is a 
hole, or exchange places with a negative particle with rate 7. At the right edge, they 
leave the system with rate f3. Negative particles enter the system from the right and 
exit at the left, in a CP-invariant fashion. In detail. 
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To implement this dynamics, we maintain a list of possible moves, draw randomly 
from that list and update according to the rates in (^. If the length of the list is I, 
continuous time, measured in units of MC steps (MCS), is incremented by 1// whenever 
a move is attempted. Taking the average of an observable then requires a weighting by 
l/l. The performance of the code can be greatly improved if we note that the particle 
density changes only at the boundaries, (III)(a)-(d) [THj. 

All simulations are performed at a = 1. Starting from an initially randomly half- 
filled system, we discard at least 10^ MCS to avoid transients, and then gather data for 
about 4 X 10^ MCS, unless otherwise noted. Smaller values of 7 require longer runs, 
up to 2.5 X 10^° MCS, to gain a clear picture. The random number generator is the 
"Mersenne twister" pTB]. 

We generate particle density histograms, Pi(p_,p+), by continuously monitoring 
the instantaneous densities of positive and negative particles, p± = (2L)~^ [\si\ ± Si], 
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in steady state. The associated configurational averages are denoted by (p+) and (p-). 

In steady state, the probabihty of a particular configuration of a finite system is 
rigorously symmetric under CP, so that Pl{p-,p+) = Pl{p+,P-) and = We 
exploit this symmetry in order to reduce noise in our statistical data, by considering the 
"folded" histogram Pl(p_, p+) + Pl(p+, In certain sectors of parameter space, CP- 
symmetry is spontaneously broken JT], in the sense that larger systems are "trapped" 
for times oc exp (const x L) in particular regions of phase space with asymmetric particle 
densities and currents. 

Even for short switching times, the particle density histograms P/^(p_,p_|_) carry 
by far the clearest signature of the transitions fT^ IT!^ . A single peak on the diagonal 
indicates a symmetric phase, while a double peak with two off-diagonal maxima signals 
a symmetry-broken phase. The transition between the two asymmetric phases is marked 
by histograms with two long ridges, one running close and nearly parallel to the p_- 
axis, and the other close to the p_|_-axis. At the transition, the global maximum suddenly 
shifts from the near (Id/ld) to the far (hd/ld) end of each ridge, like water sloshing in 
a narrow channel. Going beyond earlier studies, we quantify the histogram features 
by measuring the observables p^in and p^^j, denoting the density of the minority (pmin) 
and majority (Pmaj) species at the global maximum of Pl{p^, p+). The statistical error 
of the peak position can be determined from its scatter in independent runs. This is 
computationally very costly, so that the error bars shown for L = 500 in figure Q are 
extrapolated from a few runs at representative parameter values. They are typically of 
the size of the symbols. 

3. Results 

We first consider the L-dependence of the phase diagram, for 7 = 1. Figure C] shows 
p„ji„ and p^aj for different system sizes L = 500,2000,8000,32000 at 7 = 1. Starting 
at low P, each system size displays two branches, the upper (lower) being associated 
with p„,^j (Pmin)- The (hd/ld)- (Id/ld) transition, associated with the abrupt drop in p^^j, 
occurs at a nearly L-independent value of jSi ~ 0.2735. Even though the difference 
Pmaj — Pmin IS sigmficautly reduced at the transition, it remains nonzero above 
indicating asymmetric phases on both sides. Remarkably, only p,„^j signals the transition 
clearly while p^in is quite smooth. As we increase P further, we encounter the (Id/ld)- 
(Id-S) transition where p„,^j — p„,i„ vanishes, up to statistical fiuctuations. Labelling 
the transition point P2{L), it is obvious that P2 depends strongly on L. Our data 
for the smallest system size (500) confirm the results of [12], i.e., Pi ^ 0.274 and 
^2(500) ~ 0.284. However, the rest of our data show clearly that /32(L) is monotonically 
decreasing with L. Although the data for L = 32000 is quite noisy (with runs of 
only 2.2 x 10^ MCS overall), one can still roughly locate the second transition. In 
particular, we expect P2{L) to be just below the smallest /3 for which p^^j = p,,,;,, 
(indicated by arrows in figure Q). From these estimates of P2{L), we found that the data 
are remarkably consistent with the power law P2{L) — /?i oc The dotted line in 
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Figure 1. The position of the global maximum in the histogram Pl{p-, p+) as a 
function of /3 at 7 = 1 for different system sizes L. The filled (open) symbols show the 
density of the majority (minority) species, Pmaj and Pmin, respectively. The smallest 
value of P in the symmetric phase is marked by an arrow, indicating the onset of the 
{ld-S)-{ld/ld) transition. It moves significantly with L. In contrast, the {ld/ld)-{hd/ld) 
transition, marked by a dashed line through all four panels (smallest /3 in (Id/ld)), stays 
fixed. The dotted line (right ordinate) is a power law fit of the position of the (Id-S)- 
(Id/ld) transition, L vs. j32{L), the crosses marking the suggested transition point for 
each system size shown. 



figure^shows such a fit. Barring further surprises at even larger values of L, we take this 
as a strong indication that the (W/W)-phase does not persist in the L 00 limit. As a 
further test, we investigated the maximum of Pl{p-, p+), maXp^p^{PL(p_, p+)|p_— p+}, 
along diagonals defined by constant p_ — p+ [13 . Examples for these "ridge profiles" 
or "silhouettes" are shown in figure El For lower values of L and fixed (3 with, say, 
j3i < P < ^2(500), the ratio of peak to saddle height of this function increases with L, 
which appears to support the presence of the {ld/ld)-phase (see figure 10 in [12] )• Yet 
for sufficiently large L, this ratio decreases again and eventually drops to 1 as soon as 
P2, the position of the {ld/ld)-{ld-S) transition, has moved below the {3 considered. We 
expect to see similar behaviour for any value of /3 between 0.2735 and 0.284. 

These data provide a clear picture of the 7 = 1 model. We now turn to an 
expanded parameter space by considering smaller values of 7, reflecting congestions 
on the bridge. Based on extensive simulations (up to 5 x 10^ and 2 x 10^° MCS for the 
transient and statistics, respectively), the transitions can be investigated for 7 as low 
as 0.4. Figure El shows, similarly to figure ^ the behaviour of minority and majority 
species. Our key observation is that the width of the (/(i//(i)-phase is even further 
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Figure 2. As figure ^ the position of the global maximum in the histogram 
Pl{p-, P+) as a function of /3 for different system sizes L, but now for 7 = 0.4. Filled 
(open) symbols denote the majority (minority) density, p^^j (pmin)- The position of 
the transition {ld-S)-{ld/ld), marked by arrows, moves towards the {ld/ld)-{hd/ld) 
transition. The latter moves slightly towards larger (3 with L, the dashed line marking 
the largest (3 in [Id/ld) for the largest L. 



reduced. For 7 = 0.4 and L = 500, it is already confined to 0.175 < (3 < 0.180, 
shrinking further to 0.1763 < (3 < 0.1773 for L = 8000, with errors of at most ±0.0005. 
Unless very finely grained steps in f3 are simulated, it is easy to miss this phase entirely. 
It is also apparent that both transitions move to smaller values of /3 as 7 decreases; 
however, j32 decreases faster so that the width of the (/(i/W)-phase shrinks as 7 is 
lowered. Below 7 = 0.4 the data becomes very noisy so that only moderately large 
systems can be considered. Here, the histograms develop only short ridges, very close 
to the diagonal p+ = p_. For small /3, the global maximum is at the far end of the 
ridge, corresponding to the {hd/ld)-ph.a.se. As P increases, the ridge shortens and moves 
closer to the diagonal. During this reshaping, it becomes virtually impossible to resolve 
whether the global maximum actually moves to the near end of the ridge - which would 
be the signature of the {ld/ld)-ph.ase - or whether it just immediately appears on the 
diagonal, signalling the {ld-S)-phase. Also, the low-/3 phase is no longer characterised 
by having p,„aj > 1/2; instead, both p^^j and p^^^ are below 1/2; in this sense, the 
terminology {hd/ld) is somewhat of a misnomer. To conclude, the lowest 7 for which 
we still identify the {ld/ld)-phase with reasonable certainty is 7 = 0.4 for L = 500. The 
phase diagram in figure El summarises our findings, showing the position of the (hd/ld)- 
(Id/ld) and [ld/ld)-{ld-S) transitions in the /3,7-plane for L = 500. Below 7 = 0.4, 
we draw only a single phase boundary, separating a symmetry-broken phase from a 
symmetric one. Remarkably, the whole length of the {hd/ld)-ph.a.se boundary can be 
well approximated by the phenomenological formula 1 — a/1 — (3/ (3i with (3i = 0.2735. 
If a theoretical framework to predict this formula were developed, it would undoubtedly 
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Figure 3. Phase diagram in f3, 7 space for the bridge model at L = 500. The 
closed symbols show the transition {hd/ld)-{ld/ld), the open symbols the transition 
{ld-S)-{ld/ld), which is expected to merge with the former. The error bars indicate the 
estimated range of the transition. The dashed line is a rough fit of the first transition 
to 1 — ^Jl^^P/fii with fii = 0.2735. The two insets show typical "ridge profiles", i.e., 
the maximum of Pl{p^, p+) at given p- — p+, for the two major phases at /3 = 0.271 
and P = 0.287, both for 7 = 1 and L = 500. 



provide significant insight into the complex behaviour of this simple model. 
4. Conclusions 

Using extensive simulations for much larger system sizes than previously studied, we 
performed a finite-size analysis for a simple two-species asymmetric exclusion process 
with open boundaries. The model is reminiscent of traffic crossing a narrow bridge 
in both directions. Extending earlier work, we considered the effect of having a small 
passing probability, 7, which is the rate for two particles to exchange places. In smaller 
systems with 7 < 1, we confirm the existence of four phases. However, as L increases or 
7 decreases, the {ld/ld)-sectoT shrinks until it can no longer be resolved for 7 < 0.4 and 
L > 500. This suggests strongly that the {ld/ld)-ph.ase is a finite-size effect and that 
there is only one symmetry-breaking transition in this model. While mean-field theory 
predicts otherwise, it is often unreliable in low dimensions. Our investigation highlights 
again the importance of using very large systems in numerical studies of driven diffusive 
systems. Moreover, it suggests that scanning a wider range of control parameters may 
already reveal the same subtleties in smaller systems. Until a reliable analytic approach 
for this model is found, we believe that our study provides the best insight so far into 
the issues raised in earlier works ^21 1^ • 
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